Fluidization of wet granular matter 
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We have studied the effect of small amounts of added liquid on the dynamic behavior of a granular 
system consisting of spherical glass beads. The critical acceleration Y for fluidization is found to 
increase strongly when liquid is added. It furthermore depends on the bead radius, the height of 
the sample, and the shaking frequency. As our main result, we find that these dependencies can be 
presented in a factorized form. We present a simple model which accounts for most of our findings. 

PACS numbers: 05.65.+b,45.70.-n,45.70.Mg 



As it is generally known from everyday experience, the 
mechanical properties of a granular material change dra- 
matically if some liquid is added. The main reason is the 
internal cohesion due to capillary forces arising from liq- 
uid bridges between the grains P, 0, 0, IH H, ■ While re- 
cent years have seen considerable progress in understand- 
ing the dynamics of dry granular materials H, 0, E| i 
the physical mechanisms underlying the properties of wet 
systems still pose many unanswered questions. 

In this article we cosider granular matter in a container 
vibrated vertically with a certain frequency v and an am- 
plitude A 0. If this agitation is weak, the particles of 
the granulate remain fixed and do not move against each 
other. When the strength of the agitation is increased be- 
yond the so-called fluidization threshold by enhancing A 
or v, the granulate becomes fluidized 0. In this state the 
cohesiveness of the granulate is overcome. This enables 
a liquid-like flow of the granular matter, e.g. through 
pipe systems [llj. Furthermore there starts in general a 
convective flow resembling to that of a fluid. The form 
of the vibration used in most experiments was a sinu- 
soidal horizontal or, in particular, a vertical vibration. 

In our experiments, we have studied the effect of cap- 
illary forces, which are characteristic for wet granular 
matter, on the fluidization behavior of the granulate |27| . 
We mounted a small cylindrical glass container (diame- 
ter: 2.5 cm; height: 4 cm) filled at about three quarters 
with the granular matter consisting of glass beads on an 
inductive shaker. The container was shaken vertically 
with amplitudes between l/j,m and 1mm and frequencies 
between 20 and 333-ffz. We conducted the experiment 
with several samples composed of particles with radii 
from 138/xm to 500/im, respectively. In order to prevent 
crystallization, the beads in each sample were chosen to 
be slightly polydisperse: the spread in bead size ranged 
from 10 to 20 %. We added controlled amounts of liquid 
and shook the sample by hand for some minutes in order 
to obtain a homogeneous distribution of liquid. 

For the results presented here, we used water as the 
liquid, but similar results were obtained using nonane. 
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We can therefore exclude that leaching of ions from the 
glass, and other side effects due to the specific properties 
of water 0], play a significant role. 

We measured the fluidization threshold as follows: The 
shaking frequency v was fixed at a certain value. The 
shaking amplitude A was increased until a relative mo- 
tion of the particles was visible through the glass wall. It 
turned out that this transition was rather well defined, 
since small variations of less than 5% of the amplitude 
at the critical value were decisive for the existence of the 
fluidization behavior. From the amplitude and the fre- 
quency we calculated the fluidization amplitude: 

T := 4tt 2 

9 

where g is the gravitational acceleration. 

We found that the fluidization acceleration for wet 
granular matter depends strongly on the water content 
W, but also on the bead radius R, the sample height H 
and the shaking frequency v as follows: 

T(W, H, R, v) = r (l + G{W, H, R, v)) 

with To = 1.2 being the value for the dry case. As it 
turned out, the function G(W, H, R, v) can be factorized 
and written as 

G(W, H, R, v) = f(W)h(H)p(R)q(u) (1) 

i.e. the contribution of each single parameter can be 
treated separately. This is the main result of this paper. 
The form of the dependencies is described below. 

Dependence on water content: The fluidization ampli- 
tude r was found to depend strongly on the water con- 
tent W as shown in Fig. ^ F° r small water contents the 
fluidization acceleration T was a rapidly increasing but 
smooth function of W and went from the value for the 
dry granulate (W = 0) to a plateau value. The value 
of the water content at which this plateau was reached, 
Wb, was roughly 0.1%, as is shown in Fig. Q]for differ- 
ent frequencies. T saturated then at this plateau value 
and stayed constant till water contents of the order of 
magnitude of 10%. The plateau value depended on the 
frequency v and the bead radius R, as we will discuss 
later in the text. 
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At equilibrium, we assume that all of the liquid is in 

the bridges, w = w , where c is the coordination 

p m 

number of the network of bridges (c w 6) |20| and p p is the 

packing density of the spheres p p ps 0.64. Under dynamic 
conditions the typical bridge volume would be expected 
to be somewhat smaller than the equilibrium value. Our 
own measurements [2l] suggest that w — aR 3 W, where 
a re 0.25 0. 

We assume that the starting condition for fluidization 
due to vertical shaking is that the force F exerted from 
the vibration is big enough to break the capillary bridges 
in at least one horizontal plane in the pile. Here the 
vibrational force is given directly by Newton's law: 



0.1 



F = (r - T )Mg 



(3) 



FIG. 1: Dependence of the fluidization acceleration for differ- 
ent frequencies for R = 0.138mm on the water content. The 
fitting functions are described in the text. 



We can provide an explanation for this behavior in 
terms of the capil lary forces which have a hysteretic char- 
acter [Hill El Eg. 

For this we look at first at the capillary forces between 
two glass beads with a radius R. As it can be shown 
analytically [15| , the capillary force between two perfectly 
smooth beads depends for sufficient small liquid contents 
only on the radius R, but not on W: 

F cap = 27ri?7 

with the surface tension 7. 

On the other hand it was found experimentally EU 
that for very small water contents and beads with a rough 
surface, the capillary forces are not constant but decrease 
monotonically and smoothly to zero for W — > 0: 

F cap = 2vR 1 f{w) 

with the volume of a single capillary bridge w and / being 
a function which depends on the roughness of the sphere 
surfaces 0,0. 

The function f(w) is for perfect spheres and complete 
wetting equal to one but for a rough surface, it tends to 
zero as w — > 0. The shape of f(w) has been discussed in 
detail before Q. Based on these considerations, it may 
be approximated by 
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with wo := R8 2 characterizing the roughness amplitude, 
5. For the exponent, we have p = (2 — x)/(2 + x)> where 
X is the roughness exponent of the bead surface. For 
(typical) scratch and dig roughness as for our beads, \ is 
close to zero, and p ps 1. This corresponds to the regime 
found by Hornbaker et al. 



with the mass of the granular sample M and To := 
T(W = 0) ps 1.2 0. The mass of the (cylindrical) pile is 
given by 

M = pgp rcp HB 

with the density of the bead material (glass) p g , the ran- 
dom close packing density for beads p rcp ps 0.64 j23|, the 
height of the cylindrical sample H and its base area B. 

The expression for the force exerted from the capillar- 
ies in the relevant horizontal plane is 



F pl = ^p 2D 2R 1 f{w) 



(4) 



with the two-dimendional 'density' p 2 D of the bridges 
which have to break. 

From equating Q and J3J, one gets 



T(W) = T 1 



cf(w) 
HR 



with 



27P2 



r- . As it can be seen in Fig. Ill this can 

be used quite well to describe the Vy-dependence of T. c 
was a fitting parameter which was also a function of the 
bead radius R and the shaking frequency v (see below). 

The function f(w) connects the bridge formation di- 
rectly to the roughness. From Wo « 0.05% (see Fig. 
we get 6 ~ 500nm, which is similar to the peak-to-peak 
roughness obtained from the inspection of the beads by 
atomic force microscopy [28|. As it can be seen in Fig. 
n the dependence of the fluidization acceleration on the 
liquid content is described quite well by this theory. 

A further prediction is the dependence of T on H, 
which we could confirm at least qualitatively in addi- 
tional control measurements at 0.25mm < H < 35mm 
with R — 138/im and W = 0.13% as is shown in Fig. [21 
We could could fit the values plotted in Fig. (2J quite well 
with 



T(W,H) = To 1 + 
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FIG. 2: Dependence of the fluidization acceleration on the 
height H of the granular pile. The fitting function is described 
in the text (H c ~ 4mm). 



The finite limiting value H c may be understood as fol- 
lows. One ruptured horizontal plane may be sufficient to 
fluidize the surrounding medium only to a certain verti- 
cal distance. Hence in samples with large H, the initial 
fluidization occurs in horizontal layers with a thickness 
H c < H such that T becomes independent of H . We were 
able observe this behavior by direct optical monitoring. 
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FIG. 3: Bead radius dependence of AF := F(W, R) - V Q for 
different frequencies (circles and squares). 
The asterisks represent the radius dependence of the slope C, 
(described in the text; see Fig. |1J. 

The solid line is proportional to BT 13 with (3 = 0.5, the dotted 
line is proportional to R~ . 

Dependence on the bead radius: As already mentioned 
the dependence of T on the bead radius R proved to be 
non-trivial. We found T{W, R) — Tq to be proportional 
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0.5, as shown in Fig. [3] for H » 3cm. 



to 

We were able to describe the data with 



T(W,R) 




(6) 



The deviation of the measured values from the predic- 
tion according to Eq. [SJgets larger for smaller and smaller 
bead radii. Hence we set Rq = 0.5mm as the value where 
and |H1 yield the same value for T 29] . Therewith we 
obtained a p2D ranging between 3 and 4, which is a rea- 
sonable result. 

A possible explanation for (3 ^ 1 is that the particles of 
a wet granulate tend more and more to form (statistical 
and temporary) clusters the smaller they are. Since these 
clusters may behave like single bigger particles 0> HH i 
T(W, R) should be especially for small R smaller than 
the fluidization acceleration expected from Eq. [S] This 
is consistent with (3 < 1, as observed. 
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FIG. 4: Frequency dependence of the fluidization acceleration 
for different bead radii for W = 0.3% where r(W) has already 
saturated. The fitting functions are described in the text. 

Frequency dependence: The fluidization amplitude 
turned out to depend also on the shaking frequency 
v. In general we observed that T was a concave and 
monotonous rising function of v. This frequency depen- 
dence was most distinctive for small beads and at high 
liquid contents, although it did not depend very strongly 
on W for the bigger R. Some sample curves of this de- 
pendence are shown in Fig. ^ 

The frequency dependence could be fitted with the 
function 

r( I /) = r (i + cin(i + i y/ i y )) 



with 
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(see Eq. ^) for the applied large frequencies with the 
fitting parameter £ (and accordingly C±), as is shown in 
Fig. 0] The value of v$ turned out to be of the order of 
1Hz. The frequency dependence was more pronounced 
for the small beads, i.e. £ was a monotonous falling func- 
tion of v. It emerged that C(i?) was roughly proportional 
to R-P (see Fig. 0. 

We will now provide a qualitative explanation for this 
behavior. The amplitudes A corresponding to the fre- 
quencies v were only about 1/xm for the higher v which 
is much smaller than the radii of the glass beads. This 
means that a translational motion of the glass beads is 
more and more difficult for higher frequencies and occurs 
only through many-particle collisions. Thus the rota- 
tional degree of freedom becomes more and more impor- 
tant for the higher frequencies which leads to a concave 
form of r(i/), since rotational motion is not impeded sig- 
nificantly by small A if the particles are perfect spheres 
|3fj| . We have already mentioned the occurrence of parti- 
cle clusters in the wet granulate. Since these clusters are 
irregular, they cannot easily rotate. Hence one expects 
a more pronounced frequency dependence, i.e. a big- 
ger C, for the smaller beads where the assumed clusters 
are more pronounced (see above). This is in qualitative 
agreement with our measurements. 

The rotational motion of non-sphcrical particles in a 
granulate consisting of them is much more difficult (as for 



the case of clusters described above) and depends clearly 
on A. Hence we conducted also some experiments with 
cubes 31] and other particles like cylinders, river sand 
and foraminiferes in order to check our assumption. For 
the case of cubic particles [33] we observed indeed a very 
different behavior: The increase of T(u) was nearly lin- 
ear and quite strong |33| . The ^-dependence of T was 
for the other non-spherical particles also linear but not 
as pronounced as for the cubical case. The dependence 

Ci := in the case of cylinders was very close to 

that for the glass beads, while the cases of river sand 
and foraminiferes lied in-between the observed extremes 
given by the spheres and the cubes. Concluding it can be 
said that (i increased with the 'roughness' (or 'edginess') 
of the particles. Furthermore Ci was a monotonous rising 
function of W which is also consistent with our explana- 
tion of the ^-dependence |2J] . 

It is interesting to note that for the beads, with a 
roughness exponent x ~ 0, we get T — Tq oc ln(^), while 
for edgy grains, with \ ~ 1 we find T — To cx v . Whether 
this finding represents a coincidence or bears some fun- 
damental physics remains to be investigated. 
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